Detection of spin injection into a double quantum dot: 
Violation of magnetic-field-inversion symmetry of nuclear polarization instabilities 
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In mesoscopic systems with spin-orbit coupling, spin-injection into quantum dots at zero magnetic 
field is expected under a wide range of conditions. However, up to now, a viable approach for 
experimentally identifying such injection has been lacking. We show that electron spin injection 
into a spin-blockaded double quantum dot is dramatically manifested in the breaking of magnetic- 
field-inversion symmetry of nuclear polarization instabilities. Over a wide range of parameters, the 
asymmetry between positive and negative instability fields is extremely sensitive to the injected 
electron spin polarization and allows for the detection of even very weak spin injection. This 
phenomenon may be used to investigate the mechanisms of spin transport, and may hold implications 
for spin-based information processing. 



Time reversal symmetry is a fundamental law of na- 
ture, which places strong constraints on the types of be- 
haviors which can occur in physical systems. Based on 
the Onsager concept of microscopic reversibility close to 
equilibrium, general transport coefficients exhibit partic- 
ular symmetries under the reversal of direction of an ap- 
plied magnetic field B [T]. Such magnetic- field- inversion 
symmetry is robustly observed to high precision through- 
out a wide variety of experiments. Therefore it is quite 
remarkable to find examples of phenomena where this 
symmetry is violated. Moreover, such asymmetries can 
provide information about deviations from equilibrium. 

Over the past several years, excitement about the 
prospect of spin-based information processing has led 
many authors to consider a variety of mechanisms for 
injecting and manipulating electron spins in nanoscale 
devices. Such works have shown that, through the spin- 
orbit interaction, significant spin injection can be pro- 
duced even in the absence of an applied magnetic field [2j- 
[7]. The only restrictions on spin injection are imposed by 
the action of time reversal symmetry together with uni- 
tarity, and they allow spin injection in systems with more 
than one outgoing channel [21 El [8]. However, as shown 
in Ref. 6 , coupling to an environment which breaks uni- 
tarity allows spin injection even into a single outgoing 
channel. Experimentally, spin injection through quan- 
tum point contacts has also been reported [HI HB1 • Thus 
spin injection appears to be a generic phenomenon. 

Although less extensively studied, analogous mecha- 
nisms should lead to spin injection of electrons into quan- 
tum dots in systems with spin-orbit coupling, at B = 0. 
For quantum dots coupled to source and drain electrodes, 
unitarity is broken due to the coupling to phonons and 
Fermi reservoirs, and there is no fundamental reason to 
expect a vanishing spin injection probability. 

Although such spin injection is expected under a wide 
range of circumstances, until now a viable method for its 
detection has been missing. In this paper, we demon- 
strate that spin injection can be manifested in a dra- 
matic violation of magnetic-field-inversion symmetry in 



dc transport through spin-blockaded double quantum 
dots. Experiments in this regime |llH13j have demon- 
strated a variety of interesting nonlinear phenomena such 
as bistabilities and hysteresis, which are associated with 
the coupled dynamics of electron and nuclear spins. In 
particular, B-inversion asymmetry has apparently been 
recently observed by the Delft group [2j. To date, the 
theoretical treatment of instabilities [T5J [IB] was based 
on the assumption of completely unpolarized injected 
electron spins. Here we use an extended version of the 
model of Ref. [15] to show that spin injection breaks the 
magnetic-field-inversion symmetry of the dynamical in- 
stabilities. Furthermore, in the regime where hyperfine- 
and non-hyperfine-mediated decay rates are comparable, 
the degree of asymmetry is an extremely sensitive func- 
tion of the injected electron spin polarization. 

Why is magnetic-field-inversion symmetry violated for 
this system? The direction of the dc current flowing 
through the double dot breaks the time-reversal symme- 
try, even for very weak currents. Through spin-orbit cou- 
pling, this violation of time-reversal by the direction of 
the current is converted into spin injection into the dot. 
While the observation of spin injection per se does not re- 
veal its mechanism, by varying electrostatic gates which 
control the transport of electrons in the lead and/or bar- 
rier regions, the phenomenon can be used to investigate 
the nature of spin transport in the system. 

Note that, in the presence of a large Zeeman split- 
ting, electron spins may be injected with a high de- 
gree of polarization[T71[TE]. However, such magnetically- 
induced spin injection is symmetric in the magnetic 
field B, and therefore does not lead to magnetic-field- 
inversion asymmetry. Here we focus on the field- 
independent part of spin-injection which persists down 
to zero magnetic field. For simplicity, we further assume 
that magnetically-induced spin-injection is weak over the 
range of relevant fields. To illustrate the spin- injection in- 
duced magnetic-field-inversion-symmetry-breaking phe- 
nomenon most clearly, we start from the simple model 
of spin-blockaded transport proposed in Ref. [TJ5] and em- 
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FIG. 1: Energy level diagram of spin-blockaded double quan- 
tum dot. At large potential bias, hyperfine-assisted transi- 
tions between \T±) and the (1,1) singlet state \S) are ac- 
companied by nuclear spin flips which lead to dynamical nu- 
clear polarization and feedback through the Overhauser shift 
of the electron Zeeman energy Ez, see text for details. 

ployed to investigate nuclear spin polarization instabili- 
ties in Ref. [15] . Although the details of the results are 
model-dependent, we expect the phenomenon itself to 
persist more generally, when specific properties of realis- 
tic experimental setups are taken into account. 

The energy levels of the double quantum dot as a 
function of potential bias are depicted in FigjlJ,. For 
large positive bias, indicated by the dashed vertical line, 
the ground electronic state | S' ) is a two-electron spin- 
singlct, with large weight in the "(0,2)" orbital config- 
uration where both electrons occupy the right dot. In 
addition, the double dot supports "(1,1)" spin-singlet, 
IS*), and spin-triplet, |Ib±), states, in which electron 
density is nearly equally shared between the two dots. 
Because the triplet states are not directly coupled to the 
drain lead, the steady-state current in this spin-blockade 
regime is controlled by the rates of processes such as hy- 
perfine exchange with nuclear spins which break electron 
spin conservation. These processes, in turn, can lead to 
dynamical polarization of nuclear spins. 

We study nuclear polarization dynamics through a rate 
equation for the mean polarization x = (N + —N- ) / (JV + + 
N-), where N + (N-) is the population of nuclear spins 
oriented parallel (antiparallel) to the external magnetic 
field axis. For each electron that enters the dot, the prob- 
ability of its decay resulting in a nuclear spin flip from 
down to up (up to down) is determined by the probabil- 
ity P± of having injected the state | T + ) (| T_ )), and the 
relative competition between the hyperfine spin-exchange 
decay rate W± F and the nuclear-spin-independent escape 
rate W ln , which we take to be energy-independent. As a 
result, the nuclear polarization evolves according to 

. = h ( P+Wf F P.W™ \ 

X eN\W* F + W™ WJF + W 1 *)' [) 

where Iq is the total current through the device, e is the 
electron charge, and N = N + + AL. is the total number 
of nuclear spins in the system. The injected electron spin 
polarization is encoded in the factors P+ and P_. 



We focus on the regime of large detuning, depicted in 
FigjlJ), where the energy of the singlet level | S' ) is far 
below the energies of the blockaded triplet states. Here 
the rates W± F of nuclear spin flips arising from elastic 
hyperfine transitions between | T± ) and | S ) are obtained 
from Fermi's Golden Rule: 

where e± = s ± E z (see Figjljj), and we assume a 
Lorentzian lineshape of width 7 for the decaying sin- 
glet state. Here Ez = —^ e B + Ax is the effective Zee- 
man energy including the Overhauser shift Ax, where 
fx e = —g* [lb is electron magnetic moment with g* the ef- 
fective g-factor of the material (g* « —0.4 in GaAs), \ib 
is the Bohr magneton, B — |B| is the magnitude of the 
external magnetic field, and A is the hyperfine coupling 
strength. The matrix element Ai ~ A/yfN for electron- 
nuclear spin exchange is set by the typical scale of the 
random transverse hyperfine field. The factor 1(1 ^ 1) 
counts the available phase space for finding a properly 
oriented (down or up) nuclear spin to flip. 

We now seek the steady state values of nuclear po- 
larization, obtained by setting x — in Eq. ([!]). Trans- 
forming to the set of dimensionless parameters e = e/A, 
B = n e B/A, 7 = j/A, m = M 2 /(hW in A), and the spin- 
injection coefficient r\ = (P+ — P- ) / (P+ + P_ ) , the steady 
state values of the nuclear polarization are given by the 
third-order algebraic equation 

F(x) ee ax 3 + bx 2 + cx + d = 0, (3) 

with 

a = 1, b = rj(2i + 7m - 1) - 2B, 

c = e 2 + B 2 + 7 2 - 2e(l + r]B) + 2Brj, 

d = 2Be-r l (B 2 + e 2 + j 2 +jm). (4) 

Typically, A s» 100 fJ.eV, while the singlet-triplet splitting 
e and level width 7 can be on the ^eV scale or less. There- 
fore, below we take s, 7 <C 1. As mentioned above, we 
disregard Zeeman-splitting-induced spin injection, which 
would produce an effect even B. Thus we consider 77 as 
field independent. 

A cubic equation with real coefficients, such as that 
111 Eq.©, may have either one or three real solutions, 
depending on the values of the coefficients. Each such 
solution, which corresponds to a steady state of Eq.Q, 
can be stable or unstable, depending on whether the 
flow x tends to restore or amplify small deviations from 
the steady state. In parameter regimes where Eq.Q 
possesses two stable fixed points x = 0, the system is 
bistable and will typically exhibit hysteresis and/or pos- 
sible switching|20 . As a parameter such as the mag- 
netic field B is varied, bistability disappears at bifurca- 
tion points, where two real roots of Eq. (|3| annihilate and 
become a complex-conjugate pair. 
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FIG. 2: Polarization fixed points of Eq.|T]), satisfying Eq.|3|, 
as a function of dimensionless magnetic field strength B for 
i = 7 = 0.01, rh = 0, in absence if spin injection, rj — 0. 

A typical pattern of fixed points for systems with r\ = 
is illustrated in the instability diagram in Figj2j where we 
plot the roots of Eq.([3]) as a function of magnetic field 
B. Solid (dotted) lines indicate stable (unstable) fixed 
points. Note that in absence of spin injection, i.e for 
rj = 0, the solutions are symmetric with respect to B- 
inversion. The system exhibits bistability over a wide 
range of magnetic field strengths, with bifurcation points 
near \B\ = 0.7 where bistability disappears. 

To investigate the pattern of instabilities in more de- 
tail, we examine the discriminant of Eq.([3]), which we 
denote by A[F(x)]. For a general polynomial, the dis- 
criminant A = IX^j^j — x j) ^ s a symmetric function 
of the polynomial's roots {xi}. Each complex-conjugate 
pair of roots contributes a factor of -1 to A. Therefore 
the bifurcation points, where two real solutions merge 
and turn into a complex conjugate pair, correspond to 
the zeros (sign-changing points) of the discriminant. 

Because the discriminant A[f(x)] is a symmetric func- 
tion of the roots of F(x), it can be expressed directly as a 
polynomial in the coefficients of F(x). For a cubic poly- 
nomial of the form ([3]), the discriminant is given bv|21) 

A[F(x)} = 18abcd ~ Ab 3 d + b 2 c 2 - 4ac 3 - 27a 2 d 2 . (5) 

Thus the problem of mapping out the bifurcations of the 
fixed points of the flow x in Eq.|l]) is reduced to the 
problem of solving for the roots of A[F(x)] in Eq.([5]), 
with a, b, c, and d taken from Eq.Q. 

Because we are primarily interested in the magnetic- 
field-inversion symmetry/asymmetry of the system, we 
focus on the -B-dependence of the discriminant A[F(x)]. 
With all other parameters fixed, the equation A[F(x)} — 
yields a fifth-order polynomial in B, whose roots de- 
termine the bifurcations of the fixed points of the flow 
Q. The full expression for A[.F(x)] is quite cumber- 
some, and we do not reproduce it here. The expan- 
sion of A in the regime e, 7, to, 77 <C 1, up to third or- 
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FIG. 3: Discriminant A[F(x)}, Eq.(|5|, vs. dimensionless 
magnetic field B for e = 0.01, 7 = 0.01, r\ = 0.05, and 
rh — 0.78 (upper panel) and the corresponding fixed point 
diagram (lower panel). The bifurcation points B± at positive 
and negative field values are indicated by dashed lines. Note 
that even with small (5%) spin injection, the instability fields 
B+ and B* differ by one order of magnitude. 

der in all parameters, reads as A « A' 2 ' + A^ 3 ', with 
A< 2 ) = -4(7 2 +e 2 )B 4 + 4e 2 B 2 and A< 3 ) = -4r/(m 7 )B 5 + 
477 (m-y + 4<7 2 + 2e 2 )B 3 - 40e(e 2 + ^ 2 )B 2 - 8i]i 2 B + 32 e 3 . 
Note that 77 first appears in A*- 3 -*, which is linear in 77. 
There, rj multiplies each odd power of B, ensuring that 
the discriminant is invariant under 77 — > —rj,B — > —B. 
This is a fundamenal symmetry of the system, and holds 
to all orders. 

As demonstrated in Fig|3j for rj 7^ 0, i.e. when the 
incident current carries spin-polarization, the zeros of A 
can be highly asymmetric in 5. To explore the degree 
of magnetic-field-inversion symmetry breaking in more 
detail, in Fig|4]we plot the instability asymmetry param- 
eter B* — !(-£?* +-B* ) as a function of the spin injection 
coefficient 77, and the parameter 777 which describes the 
competition between the hyperfine transition rates W± F 
and the nuclcar-spin-indepcndcnt decay rate W m . Here 
B*, and B*_ are the upper and lower bifurcation points, 
as indicated by the dashed lines in Fig(3j For g* < 0, 
as is typical for GaAs and InAs, two other roots of the 
fifth-order equation A = are complex, while the fifth 
root is unphysical because it corresponds to a nuclear 
spin polarization of greater than 100%. 

For weak spin injection, i.e. for small 77, the instability 
asymmetry B* grows monotonically with 77. The rate at 
which B* grows with 77 is controlled by the competition 
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FIG. 4: Dependence of instability asymmetry B* = |(-B+ + 
B"L ) on the spin-injection coefficient n and the factor m, which 
describes the competition between hyperfine and nuclear- 
spin-independent decay rates. 



between hyperfine and non-hyperfine decay channels, i.e. 
by to. As illustrated by the line cuts in the upper panel, 
the asymmetry grows very sharply with r\ when the two 
escape processes compete with comparable magnitudes 
(to = 1, solid line). Note that for very large 77, we find a 
boundary beyond which all bifurcation points disappear, 
and the system is stable, being partially polarized, for all 
values of B. 

Within the quantum dot, the levels | T + ) and | T_ ) 
are defined with respect to a quantization axis which is 
directed nearly along the external field B. However, the 
polarization axis of electron spin injection is mostly de- 
termined by the spin-orbit interaction along the path be- 
tween the source and the dot. In the discussion above, we 
have implicitly assumed that injected spins were polar- 
ized along B, in which case spin injection directly leads to 
an imbalance of the probabilities P+ and P_ to load the 
I T + ) and | T_ ) states. Suppose instead that the mag- 
netic field is oriented perpendicular to the axis of electron 
spin injection. In this case, the system will on average 
have no preference for loading either | T + ) or | T_ ) , and 
therefore we would find 77 = 0. Thus we expect that, 
within the simplest model of Zeeman-field-independcnt 
spin-injection, the spin-injection coefficient 77 should vary 
like the cosine of the angle between the external field and 
the spin-injection axis. 

In conclusion, spin-orbit coupling results in spin polar- 
ization of the electrons injected from nonmagnetic elec- 
trodes into a quantum dot even in the absence of an ex- 



ternal magnetic field B. He have shown that B-inversion 
asymmetry of the nuclear polarization instabilities can 
serve as a highly sensitive tool for detecting this polar- 
ization. 
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